Universality and its Origins at the Amorphous Solidification Transition 
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Systems undergoing an equilibrium phase transition from a liquid state to an amorphous solid 
state exhibit certain universal characteristics. Chief among these are the fraction of particles that are 
randomly localized and the scaling functions that describe the order parameter and (equivalently) 
the statistical distribution of localization lengths for these localized particles. The purpose of this 
Paper is to discuss the origins and consequences of this universality, and in doing so, three themes are 
explored. First, a replica-Landau- type approach is formulated for the universality class of systems 
that are composed of extended objects connected by permanent random constraints and undergo 
amorphous solidification at a critical density of constraints. This formulation generalizes the cases of 
randomly cross-linked and end-linked macromolecular systems, discussed previously. The universal 
replica free energy is constructed, in terms of the replica order parameter appropriate to amorphous 
solidification, the value of the order parameter is obtained in the liquid and amorphous solid states, 
and the chief universal characteristics are determined. Second, the theory is reformulated in terms 
of the distribution of local static density fluctuations rather than the replica order parameter. It is 
shown that a suitable free energy can be constructed, depending on the distribution of static density 
fluctuations, and that this formulation yields precisely the same conclusions as the replica approach. 
Third, the universal predictions of the theory are compared with the results of extensive numerical 
simulations of randomly cross-linked macromolecular systems, due to Barsky and Plischke, and 
excellent agreement is found. 
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I. INTRODUCTION 



During the last decade there has been an ongoing effort 
to obtain an ever more detailed understanding of the be- 
havior of randomly cross-linked macromolecular systems 
near the vulcanization transition 0-^. This effort has 
been built from two ingredients: (i) the Deam-Edwards 
formulation of the statistical mechanics of polymer net- 
works 1^]; and (ii) concepts and techniques employed in 
the study of spin glasses As a result, a detailed mean- 
field theory for the vulcanization transition — an example 
of an amorphous solidification transition — has emerged, 
which makes the following predictions: (i) For densities 
of cross-links smaller than a certain critical value (on 
the order of one cross-link per macromolecule) the sys- 
tem exhibits a liquid state in which all particles (in the 
context of macromolecules, monomers) are delocalized. 
(ii) At the critical cross-link density there is a continuous 
thermodynamic phase transition to an amorphous solid 
state, this state being characterized by the emergence of 
random static density fluctuations, (iii) In this state, a 
nonzero fraction of the particles have become localized 
around random positions and with random localization 
lengths (i.e., r.m.s. displacements), (iv) The fraction of 
localized particles grows linearly with the excess cross- 
link density, as does the characteristic inverse square lo- 
calization length, (v) When scaled by their mean value, 
the statistical distribution of localization lengths is uni- 
versal for all near-critical cross-link densities, the form 
of this scaled distribution being uniquely determined by 



a certain integro-differential equation. For a detailed re- 
view of these results, see Ref. Q; for an informal discus- 
sion, see Ref. Q. 

In the course of the effort to understand the vulcaniza- 
tion transition for randomly cross-linked macromolecular 
systems, it has become clear that one can also employ 
similar approaches to study randomly end-linked macro- 
molecular systems and also randomly cross- linked 
manifolds (i.e., higher dimensional objects) in each 
case, a specific model has been studied. For example, in 
the original case of randomly cross-linked macromolecu- 
lar systems, the macromolecules were modeled as flexible, 
with a short-ranged excluded- volume interaction, and the 
cross-links were imposed at random arc-length locations. 
On the other hand, in the case of end-linked systems, 
although the excluded-volume interaction remained the 
same, the macromolecules were now modeled as either 
flexible or stiff, and the random linking was restricted to 
the ends of the macromolecules. Despite the differences 
between the unlinked systems and the styles of linking, 
in all cases identical critical behavior has been obtained 
in mean-field theory, right down to the precise form of 
the statistical distribution of scaled localization lengths. 

Perhaps even more strikingly, in extensive numeri- 
cal simulations of randomly cross-linked macromolecu- 
lar systems, Barsky and Plischke [p^ , p^ have employed 
an off-lattice model of macromolecules interacting via a 
Lennard-Jones potential. Yet again, an essentially identi- 
cal picture has emerged for the transition to and proper- 
ties of the amorphous solid state, despite the substantial 
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differences between physical ingredients incorporated in 
the simulation and in the analytical theory. 

In the light of these observations, it is reasonable to ask 
whether one can find a common theoretical formulation 
of the amorphous solidification transition (of which the 
vulcanization transition is a prime example) that brings 
to the fore the emergent collective properties of all these 
systems that are model- independent, and therefore pro- 
vide useful predictions for a broad class of experimentally 
realizable systems. The purpose of this Paper is to ex- 
plain how this can be done. In fact, we approach the 
issue in two distinct (but related) ways, in terms of a 
replica order parameter and in terms of the distribution 
of random static density fluctuations, either of which can 
be invoked to characterize the emergent amorphous solid 
state. 

The outline of this Paper is as follows. In Sec. |l| we 
construct the universal replica Landau free energy of the 
amorphous solidification transition. In doing this, we 
review the replica order parameter for the amorphous 
solid state and discuss the constraints imposed on the 
replica Landau free energy by (a) symmetry considera- 
tions, (b) the smallness of the fraction of particles that 
are localized near the transition, and (c) the weakness 
of the localization near the transition. In Sec. HI we in- 
voke a physical hypothesis to solve the stationarity con- 
dition for the replica order parameter, thereby obtaining 
a mean-field theory of the transition. We exhibit the uni- 
versal properties of this solution and, in particular, the 
scaling behavior of certain central physical quantities. In 
Sec. IV we describe an alternative approach to the amor- 



phous solidification transition, in which we construct and 
analyze the Landau free energy expressed in terms of the 
distribution of static density fiuctuations. Although we 
shall invoke the replica approach in the construction of 
this Landau free energy, its ultimate form does not refer 
to replicas. As we show, however, the physical content 
of this Landau theory is identical to that of the replica 
Landau theory addressed in Sees. || and III, In Sec. 
we exhibit the predicted universality by examining the 
results of extensive numerical simulations of randomly 
cross-linked macr om olecular networks, due to Barsky and 
Plischke. In Sec. VI we give some concluding remarks. 



II. UNIVERSAL REPLICA FREE ENERGY FOR 
THE AMORPHOUS SOLIDIFICATION 
TRANSITION 

We are concerned, then, with systems of extended ob- 
jects, such as macromolecules, that undergo a transi- 
tion to a state characterized by the presence of random 
static fluctuations in the particle-density when subjected 
to a sufficient density of permanent random constraints 
(the character and statistics of which constraints preserve 
translational and rotational invariance). In such states, 
translational and rotational symmetry are spontaneously 
broken, but in a way that is hidden at the macroscopic 



level. We focus on the long wavelength physics in the 
vicinity of this transition. 

In the spirit of the standard Landau approach, we en- 
visage that the replica technique has been invoked to 
incorporate the consequences of the permanent random 
constraints, and propose a phenomenological mean-field 
replica free energy, the n —>^ limit of which gives the 
disorder-averaged free energy, in the form of a power se- 
ries in the replica order parameter. We invoke symmetry 
arguments, along with the recognition that near to the 
transition the fraction of particles that are localized is 
small and their localization is weak. The control pa- 
rameter e is proportional to the amount by which the 
constraint density exceeds its value at the transition. As 
we shall see, the stationarity condition for this general, 
symmetry-inspired Landau free energy is satisfied by pre- 
cisely the order-parameter hypothesis that exactly solves 
the stationarity conditions derived from semi-microscopic 
models of cross-linked and end-linked macromolecules. 
From the properties of this solution we recover the pri- 
mary features of the liquid-amorphous solid transition. 

In a system characterized by static random density- 
fiuctuations, one might naively be inclined to use the 
particle-density as the order parameter. However, the 
disorder-averaged particle density cannot detect the tran- 
sition between the liquid and the amorphous solid states, 
because it is uniform (and has the same value) in 
both states: a subtler order parameter is needed. As 
shown earlier, for the specific cases of randomly cross- 
linked and end-linked macromolecular systems, 
the appropriate order parameter is instead: 
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where N is the total number of particles, (with i — 
1, . . . , N) is the position of particle i, the wave- vectors 
k^, k^, • • • , kf are arbitrary, {■ ■ ■)x denotes a thermal av- 
erage for a particular realization x of the disorder, and 
[• • •] represents averaging over the disorder. 

We make the Deam-Edwards assumption [|j that the 
statistics of the disorder is determined by the instanta- 
neous correlations of the unconstrained system. (It is 
as if one took a snapshot of the system and, with some 
nonzero probability, added constraints only at those lo- 
cations where two particles are in near contact.) This 
assumption leads to the need to work with the n — > 
limit of systems of rt + 1, as opposed to n, replicas. The 
additional replica, labeled by a = 0, represents the de- 
grees of freedom of the original system before adding the 
constraints, or, equivalently, describes the constraint dis- 
tribution. 

We assume, for the most part, that the free energy is 
invariant under the group Sn+i of permutations of the 
n + 1 replicas. In this replica formalism, the replica or- 
der parameter turns out to be 
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Here, hatted vectors denote replicated collections of vec- 
tors, viz., V = {v°, v-^, • • • , v"-}, their scalar product be- 
ing V ■ w = J2a=o ■ ^'^^ (' ' ')n+i denotes an av- 
erage for an effective pure (i.e., disorder- free) system of 
n + 1 coupled replicas of the original system. We use the 
terms one-replica sector and higher-replica sector to refer 
to replicated vectors with, respectively, exactly one and 
more than one replica a for which the corresponding vec- 
tor k" is nonzero. In particular, the order parameter in 
the one-replica sector reduces to the disorder-averaged 
mean particle-density, and plays only a minor role in 
what follows. The appearance of n -I- 1 replicas in the 
order parameter allows one to probe the correlations be- 
tween the density fluctuations in the constrained system 
and the density fluctuations in the unconstrained one. 

We first study the transformation properties of the 
order parameter under translations and rotations, and 
then make use of the resulting information to determine 
the possible terms appearing in the replica free energy. 
Under independent translations of all the replicas, i.e., 
cf — » of -I- a", the replica order parameter, Eq. (||), trans- 
forms as 
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(3) 



Under independent rotations of the replicas, defined by 
Rv = {i?"v°, • • ■ , i?"v"} and cf R°'cf, the order pa- 
rameter transforms as 



^k 



^'k 
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(4) 



As discussed in detail in [Q , because we are concerned 
with the transition between liquid and amorphous solid 
states, both of which have uniform macroscopic density, 
the one-replica sector order parameter is zero on both 
sides of the transition. This means that the sought free 
energy can be expressed in terms of contributions refer- 
ring to the higher-replica sector order parameter, alone. 

We express the free energy as an expansion in (inte- 
gral) powers of the replica order parameter fi^, , retaining 
the two lowest possible powers of 57^, which in this case 
are the square and the cube. We consider the case in 
which no external potential couples to the order parame- 
ter. Hence, there is no term linear in the order parameter. 
We make explicit use of translational symmetry, Eq. (||), 
and thus obtain the following expression for the replica 
free energy (per particle per space dimension) JF„({51^}) 



nd.F„({%}) =^^.g2(fc)|f^^P 

Here, the symbol ^ denotes a sum over replicated vectors 
k lying in the higher-replica sector, and we have made ex- 
plicit the fact that the right hand side is linear in n (in the 



n — > limit) by factoring n on the left hand side. In a mi- 
croscopic approach, the functions 52 (A:) and (?3(fci, k2, k^) 
would be obtained in terms of the control parameter e, 
together with density correlators of an uncross-linked liq- 
uid having interactions renormalized by the cross-linking. 
Here, however, we will ignore the microscopic origins of 
(72 and (73, and instead use symmetry considerations and 
a long-wavelength expansion to determine only their gen- 
eral forms. In the saddle-point approximation, then, the 
disorder-averaged free energy / (per particle and space 
dimension) is given by 0: 



f = lim min J-„(i fir}) ■ 



(6) 



Bearing in mind the physical notion that near the tran- 
sition any localization should occur only on long length- 
scales, we examine the long wavelength limit by also 
performing a low-order gradient expansion. In the term 
quadratic in the order parameter we keep only the leading 
and next-to-leading order terms in k] in the cubic term 
in the order parameter we keep only the leading term in 
k. Thus, the function (73 in Eq. (^ is replaced by a con- 
stant and the function 32 is expanded to quadratic order 
in k. By analyticity and rotational invariance, 32 can 
only depend on {k°, . . . , k"} via {|k°p, . . . , |k"p}, and, 

in particular, terms linear in k are excluded. In addition, 
by the permutation symmetry among the replicas, each 
term |k"p must enter the expression for 32 with a com- 
mon prefactor, so that the dependence is in fact on k^. 
Thus, the replica free energy for long-wavelength density 
fluctuations has the general form: 

"•^X^feifc^fcs^^l ^fcs ^kl+k2 + k3,6 ■ C^) 

To streamline the presentation, we take advantage of the 
freedom to rescale J-'n, e and k, thus setting to unity the 
parameters a, b and c. Thus, the free energy becomes 
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By taking the first variation with respect to il ^, we 
obtain the stationarity condition for the replica order pa- 
rameter: 



= nd 



5Tn 



(9) 



kik2 



This self-consistency condition applies for all values of k 
lying in the higher-replica sector. 
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III. UNIVERSAL PROPERTIES OF THE ORDER 
PARAMETER IN THE AMORPHOUS SOLID 
STATE 

Generalizing what was done for cross-linked and end- 
linked macromolecular systems, we hypothesize that the 
particles have a probability q of being localized (also 
called the "gel fraction" in the context of vulcanization) 
and 1 — q oi being delocalized, and that the localized 
particles are characterized by a probability distribution 
2C~'^p(C~^) for their localization lengths ^. Such a char- 
acterization weaves in the physical notion that amor- 
phous systems should show a spectrum of possibilities 
for the behavior of their constituents, and adopts the 
perspective that it is this spectrum that one should aim 
to calculate. This hypothesis translates into the following 
expression for the replica order parameter 



(10) 



where we have used the notation k = X]q=o ■ '^^'^ first 
term on the right hand side term represents delocalized 
particles, and is invariant under independent translations 
of each replica (cf. Eq. In more physical terms, this 
corresponds to the fact that not only the average particle 
density but the individual particle densities are transla- 
tionally invariant. The second term represents particles 
that are localized, and is only invariant under common 
translations of the replicas (i.e., translations in which 
a" = a for all a). This corresponds to the fact that 
the individual particle density for localized particles is 
not translationally invariant, so that translational invari- 
ance is broken microscopically, but the average density 
remains translationally invariant, i.e. the system still is 
macroscopically translationally invariant (MTI) . 

By inserting the hypothesis ( p^ into the stationarity 
condition (^, and taking the n ^ limit, we obtain 



= Sf'^ 

k. 
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In the limit fc^ — )■ 0, the equation reduces to a condition 
for the localized fraction q, viz.. 







-2qe + 3q^ 



(12) 



For negative or zero e, corresponding to a constraint den- 
sity less than or equal to its critical value, the only phys- 
ical solution is g = 0, corresponding to the liquid state. 
In this state, all particles are delocalized. For positive 
e, corresponding to a constraint density in excess of the 
critical value, there are two solutions. One, unstable, 
is the continuation of the liquid state q = 0; the other, 
stable, corresponds to a nonzero fraction, 



(13) 



being localized. We identify this second state as the 
amorphous solid state. From the dependence of the local- 
ized fraction q on the control parameter e and the form 
of the order parameter Eq. (^0|) we conclude that there 
is a continuous phase transition between the liquid and 
the amorphous solid states at e = 0, with localized frac- 
tion exponent (3=1 (i.e., the classical exponent 
It is worth mentioning that microscopic approaches go 
beyond this linear behavior near the transition, yielding 
a transcendental equation for q{e), valid for all values of 
the control parameter e; see Ref. [HJJ]. From Eq. || it is 
evident that the liquid state is locally stable (unstable) 
for negative (positive) e: the eigenvalues of the resulting 
quadratic form are given by X{k) = — e + k'^ /2. 

Now concentrating on the amorphous solid state, by 
inserting the value of the localized fraction into Eq. ([ll|), 
we obtain the following integro-differential equation for 
the probability distribution for the localization lengths: 



dTlp{Tl)p{T - n). (14) 



All parameters can be seen to play an elementary role 
in this equation by expressing p(t) in terms of a scaling 
function: 

p(r) = (2/e)7r(0); r={e/2)9. (15) 
Thus, the universal function 7t{9) satisfies 



(16) 



Solving this equation, together with the normalization 
condition 1 = d9 ^{B), one finds the scaling function 
shown in Rcfs. [^|J^]. The function 7r(6') has a peak at 
~ 1 of width of order unity, and decays rapidly both 
as — >• and 6 ^ 00. By combining these features of 
t:{9) with the scaling transformation ( p^ ) we conclude 
that the typical localization length scales as e~^/^ near 
the transition. The order parameter also has a scaling 
form near the transition: 



% = (1 - (2^/3))'5fe,o + (2e/3)s(^luj[^2kye), 
Lu{k) = / d6l7r(0)e-'='/2^ (17) 



Equation ( |16[ ) and the normalization condition on 7r(0) 
are precisely the conditions that determine the scaling 
function for the cross-linked and end-linked cases; they 
are discussed extensively, together with the properties of 
the resulting distribution of localization lengths and or- 
der parameter, in Refs. 

As discussed in this section, the localized fraction g(e) 
and the scaled distribution of inverse square localization 
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lengths 7r(0) are universal near the transition. We now 
discuss this issue in more detail. 

First, let us focus at the mean- field level. Recall the 
mean-field theory of ferromagnetism |T5| and, in partic- 
ular, the exponent which characterizes the vanishing 
of the magnetization density order parameter (from the 
ferromagnetic state) as a function of the temperature at 
zero applied magnetic field. The exponent /3 takes on the 
value of 1/2, regardless of the details of the mean- field 
theory used to compute it. The functions q{e), T^iO) and 
uj{k) are universal in the same sense. The case of g(e) 
is on essentially the same, standard, footing as that of 
the magnetization density. What is not standard, how- 
ever, is that describing the (equilibrium) order parame- 
ter is a universal scaling function, Lu{k) [or, equivalently, 
7r(0)] that is not a simple power law. This feature arises 
because the usual presence of fields carrying internal in- 
dices, such as cartesian vector indices in the case of fer- 
romagnetism, is replaced here by the external continuous 
replicated wave vector variable k. There are two facets 
to this universal scaling behavior of the order parame- 
ter. First, for systems differing in their microscopic de- 
tails and their constraint densities there is the possibility 
of collapsing the distribution of localization lengths on 
to a single function, solely by rescaling the independent 
variable. Second, there is a definite prediction for the 
dependence of this rescaling on the control parameter e. 

Now, going beyond the mean-field level, in the con- 
text of vulcanization de Gennes ||l^ has shown that the 
width of the critical region, in which fluctuations domi- 
nate and mean-field theory fails, vanishes in the limit of 
very long macromolecules in space-dimension d = 3 or 
higher. Thus, one may anticipate that for extended ob- 
jects the mean-field theory discussed here will be valid, 
except in an exceedingly narrow region around the tran- 
sition. Nevertheless, if — as is usually the case — the effec- 
tive hamiltonian governing the fluctuations is the Lan- 
dau free energy then the universality discussed here is 
expected to extend, mutatis mutandis, into the critical 
regime. 

That the amorphous solid state given by Eq. (|l^) is 
stable with respect to small perturbations (i.e., is locally 
stable) can be shown by detailed analysis. Moreover, as 
we shall see in Sec. it yields predictions that are in 
excellent agreement with subsequent computer simula- 
tions. However, there is, in principle, no guarantee that 
this state is globally stable (i.e., that no states with lower 
free energy exist). 

Up to this point we have assumed that the free energy 
is invariant under interchange of all replicas, including 
the one representing the constraint distribution (a = 0), 
with any of the remaining n, i.e. that the free energy is 
symmetric under the group Sn+i of permutations of all 
n + 1 replicas. This need not be the case, in general, as 
the system can be changed, (e.g., by changing the tem- 
perature) after the constraints have been imposed. In 
this latter case, the free energy retains the usual Sn sym- 



metry under permutations of replicas a = 1, . . . ,n. The 
argument we have developed can be reproduced for this 
more general case with only a minor change: in the free 
energy, we can no longer invoke Sn+i symmetry to argue 
that all of the |k"p must enter the expression for 32 with 
a common prefactor. Instead, we only have permutation 
symmetry among replicas a = l,...,n and, therefore, 
the prefactor b for all of these replicas must be the same, 
but now the prefactor bo for replica a = can be differ- 
ent. This amounts to making the replacement 



bob- 



En 
k" 
a=l ' ' 



(18) 



in the free energy. Both the rest of the derivation and 
the results are unchanged, except that k^ needs to be re- 
placed by fc^, throughout. We mention, in passing, that 
no saddle points exhibiting the spontaneous breaking of 
replica permutation symmetry have been found, to date, 
either for systems with Sn+i or Sn symmetry. 



IV. FREE ENERGY IN TERMS OF THE 
DISTRIBUTION OF STATIC DENSITY 
FLUCTUATIONS 

The aim of this section is to construct an expression for 
the disorder-averaged Landau free energy for the amor- 
phous solidification transition, J-, in terms of the dis- 
tribution of static density fluctuations. We present this 
approach as an alternative to the strategy of construct- 
ing a replica free energy J^n in terms of the replica order 
parameter f2. In the familiar way, the equilibrium state 
will be determined via a variational principle: SJ- = 
and S^T > 0. What may be less familiar, however, is 
that in the present setting the independent variables for 
the variation (i.e. the distribution of static density fluc- 
tuations) themselves constitute a functional. 

Our aim, then, is to work not with the replica order pa- 
rameter flj^ , but instead with the disorder-averaged prob- 
ability density functional for the random static density 
fluctuations ]l6|,^, A/'({/9k}), which is defined via 



Af({a,)) 



i=l k 



Pk 



(exp (ik • Ci)), 



(19) 



Here, Yik denotes the product over all d-vectors k, and 
the Dirac (5-function of complex argument Sc{z) is defined 
by dc{z) = S(Re z) S{Im. z), where He z and Im^; respec- 
tively denote the real and imaginary parts of the complex 
number z. From the definition of A/'({pk}), we see that 
p_k = pj^ and po = 1- Thus we can take as independent 
variables pk foi' all c?-vectors k in the half-space given by 
the condition k • n > for a suitable unit d-vector n. In 
addition, A/'({pk}) obeys the normalization condition 



I I?pAA({pk}) = 1. 



(20) 
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It is straightforward to check that, for any particular pos- 
itive integer g, the rephca order parameter flf, is a g^^ 
moment of M{{pk}) ■ 



Vp 7V({pk}) Pkl Pk2 • • • PkS = f^kl,k2,...,k9, (21) 



We use T>p to denote the measure Yi^dKe pit^dlm p^. 

The merit of the distribution functional A/'({pk}) over 
the replica order parameter Qj^ is that, as we shall soon 
see, it allows us to formulate a Landau free energy for 
the amorphous solidification transition, depending on 
A/'({pk}), hi which replicated quantities do not appear, 
while maintaining the physical content of the theory. At 
the present time, this approach is not truly independent 
of the replica approach, in the following sense: we employ 
the replica approach to derive the free energy, Eq. (0), 



and only then do we transform from the language of order 
parameters to the language of the distribution of static 
density fluctuations. We are not yet in possession of ei- 
ther an analytical scheme or a set of physical arguments 
that would allow us to construct the Landau free en- 
ergy directly. Nevertheless, we are able, by this indirect 
method, to propose a (replica-free) free energy, and also 
to hypothesize (and verify the correctness of) a stationary 
value of A/'({/9k})- It would, however, be very attractive 
to find a scheme that would allow us to eschew the replica 
approach and work with the distribution of static density 
fluctuations from the outset. 

To proceed, we take the replica Landau free energy 



Eq. (g|), in terms of the replica order parameter f^^, and 
replace flf. by its expression in terms of the {n+ 1)"^ mo- 
ment of AA({pk}). Thus, we arrive at the replica Landau 
free energy: 



ndTn ^ ( - 2 + {3 - e) J Ppi AA({pi,k}) Pp2 A/'({p2,k}) [Y.^ Pl,kP2,-k 

-t-i(n-Hl) J I?PlAA({pi,k})Pp2AA({p2,k}) (5IkfcVi.kP2,-k) (l]^Pl,kP2,-k 
- / PpiA/'({pi,k})2?P2A/'({p2,k})l?P3AA({p3,k}) , Pl,kiP2,k2P3,-ki-k2 I 

In order to obtain the desired (replica- independent) free energy we take the limit n of Eq. ([2^): 

dJ^ = d hni :r„ (3 - e) ^ Ppi A/'({pi,k}) I'p2 A/'({/92,k}) (^^^Pi,i<- P2.-kj In (^^^ pi,k P2,-k 
+ - / Vp,M{{piM})'Dp2U{{p2M}) (5Ik^''°i'kP2.-k) 

Pl,kP2,-k 



(22) 



-J PpiA/'({pi,k})2?P2A/'({p2,k})Pp3AA({p3,k}) (Ek,,k. 



Pl,ki P2,k2 P3,-ki-k 







Xing: 



'ki.k: 



Pl,ki P2,k2 P3,-ki-ks 



(23) 



In deriving the above free energy we have employed the physical fact that the average particle-density is uniform. 
In other words, the replica order parameter is zero if all but one of the replicated wave vectors is nonzero which, 
translated in the language of the distribution of static density fluctuations, means that the first moment of the static 
density distribution equals (5k, o- It is worth noting that, within this formalism, the replica limit can already be taken 
at the level of the free energy, prior to the hypothesizing of an explicit form for the stationary value of the order 
parameter. On the one hand, this is attractive, as it leads to a Landau theory in which replicas play no role. On the 
other hand, the approach is, at present, restricted to replica-symmetric states. 

We now construct the self-consistency condition that follows from the stationarity of the replica-independent free 
energy. We then proceed to solve the resulting functional equation exactly, by hypothesizing a solution having precisely 
the same physical content as the exact solution of the replica self-consistency condition discussed in Sec. |l|. 

To construct the self-consistency condition for A/'({pk}) it is useful to make two observations. First, A/'({pk}) 
obeys the normalization condition (|20|). This introduces a constraint on the variations of A/'({pk}) which is readily 
accounted for via a Lagrange's method of undetermined multipliers. Second, as mentioned above, not all the variables 
{pk} are independent: we have the relations po = 1 and p_k = Pk- I'^ principle, one could proceed by defining a 
new distribution that only depends on the independent elements of {pk}, and re-express the free energy in terms of 
this new distribution. However, for convenience we will retain Ar({pk}) as the basic quantity to be varied, and bear 
in mind the fact that po = 1 and p_k = Pk- By performing the constrained variation of dJ- with respect to the 
functional, A/'({pk}) 







^AA({pk}) 



^ + ^J ^^Pl^f{{plM}) 



(24) 
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where A is the undetermined multiplier, we obtain the self-consistency condition obeyed by A/'({/9k}): 
= Xd + 2{S-e) J VpiM{{p,M}) (Ek^kPi,-k) ln(^^PkPi,- 



(25) 



To solve this self-consistency condition for A/'({pk}) we import our experience with the replica approach, thereby 
constructing the normalized hypothesis 



/J /'OC 
- dTp{r)l[S,{p^-e^^- 



k-fc^/2T\ 



(26) 



in which q (which satisfies < g < 1) is the localized fraction and p(t) (which is regular and normalized to unity) is 
the distribution of localization lengths of localized particles. It is straightforward to show that by taking the (n -I- 1)*** 
moment of A/'({pk}) we recover the self-consistent form of the replica order parameter, Eq. (|o|). 

By inserting the hypothesis ( p6| ) into Eq. (|2^), making the replacement po — * 1, and performing some algebra, the 
self-consistency condition takes the form 







dc 



k^O 

x{2q{~e + 3q)p{T) 



'k^O 



/OkG 



-ic-k-fe^/2T 



q — {2TMr))-Sq^ 



-dq' 



dTi p{ti) p{t - Tl) 
dTi p{ti) dT2 p{t2) ln|y2/'^riT2/27re(Ti -t-r2)| -f Ad, 



in terms of the undetermined multiplier A. To determine A we insert the choice pk ~ <^k,o, which yields 

3 



A = -c 



^ dn p{ti) dT2p{T2) In {y2/'iriT2/27re(Ti + T2)} 



(27) 



(28) 



By using this result to eliminate A from the self-consistency condition, and observing that this condition must be 
satisfied for arbitrary {/Ck}, we arrive at a condition on q and p{t), viz., 



= 2q{~e + 3q)p{T) 



q—(2T'^p{T)j -3q^y dTip{Ti)p{T - Tl). 



(29) 



We integrate this equation over all values of r and use 
the normalization condition on pir) to arrive at the same 
equation relating q and e as was found in Eq. (^2|) of 
the previous section, the appropriate solution of which 
is given by g = 2e/3, i.e., Eq. ([l3|). Finally, we use this 
result for q to eliminate it from Eq. ([29|), thus arriving at 
the same self-consistency condition on p{t) as was found 
in Eq. (^4|) of the previous section. Thus, we see that 
these conditions, one for q and one for pir), are precisely 
the same as those arrived at by the replica method dis- 



cussed in Sec. Ill 



V. COMPARISON WITH NUMERICAL 
SIMULATIONS: UNIVERSALITY EXHIBITED 



The purpose of the present section is to examine the 
conclusions of the Landau theory, especially those con- 
cerning universality and scaling, in the light of the ex- 
tensive molecular dynamics simulations, performed by 
Barsky and Plischke These simulations address 

the amorphous solidification transition in the context of 
randomly cross-linked macromolecular systems, by using 
an off-lattice model of macromolecules interacting via a 
Lennard-Jones potential. 
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0.5^ 



q(n) 




FIG. 1. Localized fraction g as a function of the number 
of cross-links per macromolecule n, as computed in molecu- 
lar dynamics simulations by Barsky and Plischke (1997, un- 
published). L is the number of monomers in each macro- 
molecule; A*' is the number of macromolecules in the system. 
The straight line is a linear fit to the A'^ = 200 data. Note 
the apparent existence of a continuous phase transition near 
n = 1, as well as the apparent linear variation of q with n, 
both features being consistent with the mean-field descrip- 
tion. 



It should be emphasized that there are substantial dif- 
ferences between ingredients and calculational schemes 
used in the analytical and simulational approaches. In 
particular, the analytical approach: (i) invokes the 
replica technique; (ii) retains interparticle interactions 
only to the extent that macroscopically inhomogeneous 
states are disfavored (i.e., the one- replica sector remains 
stable at the transition); (iii) neglects order-parameter 
fluctuations, its conclusions therefore being independent 
of the space-dimension; and (iv) is solved via an Ansatz, 
which is not guaranteed to capture the optimal solution. 



o n=1.0 
o n=1.025 

n=1.05 
A n=1.075 




FIG. 2. Unsealed probability distribution of localiza- 
tion lengths ^ (in units of the linear system size), as com- 
puted in molecular dynamics simulations by Barsky and Plis- 
chke (1997, unpublished). In the simulations the number of 
segments per macromolecule is 10; and the number of macro- 
molecules is 200. 

Nevertheless, and rather strikingly, the simulations 
yield an essentially identical picture for the transition 
to and properties of the amorphous solid state, inas- 
much as they indicate that (i) there exists a (cross-link- 
density controlled) continuous phase transition from a 
liquid state to an amorphous solid state; (ii) the criti- 
cal cross-link density is very close to one cross-link per 
macromolecule; (iii) q varies linearly with the density of 
cross-links, at least in the vicinity of this transition (see 
Fig. 0); (iv) when scaled appropriately (i.e., by the mean 
localization length), the simulation data for the distribu- 
tion of localization lengths exhibit very good collapse on 
to a universal scaling curve for the several (near-critical) 
cross-link densities and macromolecule lengths consid- 
ered (see Figs. || and ||); and (v) the form of this uni- 
versal scaling curve appears to be in remarkably good 
agreement with the precise form of the analytical predic- 
tion for this distribution. 



n=1.0 
n=1.025 
n=1.05 
n=1.075 

P, 




FIG. 3. Probability distribution (symbols) Ps of localiza- 
tion lengths ^, scaled with the sample-average of the localiza- 
tion lengths ^av, as computed in molecular dynamics simula- 
tions by Barsky and Plischke (1997, unpublished). Note the 
apparent collapse of the data on to a single universal scaling 
distribution, as well as the good quantitative agreement with 
the mean-field prediction for this distribution (solid line). In 
the simulation the number of segments per macromolecule is 
10; and the number of macromolecules is 200. The mean-field 
prediction for Psc(C/Cav) is obtained from the universal scaling 
function n{9) by -Psc(y) ~ {2s/y'^) 7r(s/j/^), where the constant 
s ~ 1.224 is fixed by demanding that dyyPsc{y) ~ 1. 

It should not be surprising that by focusing on univer- 
sal quantities, one finds agreement between the analyt- 
ical and computational approaches. This indicates that 
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the proposed Landau theory does indeed contain the es- 
sential ingredients necessary to describe the amorphous 
sohdification transition. 

Let us now look more critically at the comparison be- 
tween the results of the simulation and the mean-field 
theory. With respect to the localized fraction, the Lan- 
dau theory is only capable of showing the linearity of the 
dependence, near the transition, on the excess cross-link 
density, leaving undetermined the proportionality factor. 
The simulation results are consistent with this linear de- 
pendence, giving, in addition, the amplitude. There are 
two facets to the universality of the distribution of lo- 
calization lengths, as mentioned in Sec. [II. First, that 
the distributions can, for different systems and different 
cross-link densities, be collapsed on to a universal scal- 
ing curve, is verified by the simulations, as pointed out 
above. Second, the question of how the scaling param- 
eter depends on the excess cross-link density is difficult 
to address in current simulations, because the dynamic 
range for the mean localization length accessible in them 
is limited, so that its predicted divergence at the transi- 
tion is difficult to verify. 



VI. SUMMARY AND CONCLUDING REMARKS 

To summarize, we have proposed a replica Landau 
free energy for the amorphous solidification transition. 
The theory is applicable to systems of extended objects 
undergoing thermal density fluctuations and subject to 
quenched random translationally-invariant constraints. 
The statistics of the quenched randomness are deter- 
mined by the equilibrium density fluctuations of the un- 
constrained system. We have shown that there is, generi- 
cally, a continuous phase transition between a liquid and 
an amorphous solid state, as a function of the density 
of random constraints. Both states are described by ex- 
act stationary points of this replica free energy, and are 
replica symmetric and macroscopically translationally in- 
variant. They differ, however, in that the liquid is trans- 
lationally invariant at the microscopic level, whereas the 
amorphous solid breaks this symmetry. 

We have also shown how all these results may be recov- 
ered using an alternative formulation in which we focus 
less on the replica order parameter and more on the dis- 
tribution of random static density fluctuations. In par- 
ticular, we construct a representation of the free energy 
in terms of this distribution, and solve the resulting sta- 
tionarity condition. 

Lastly, we have examined our results in the light of the 
extensive molecular dynamics simulations of randomly 
cross-linked macromolecular systems, due to Barsky and 
Plischke. Not only do these simulations support the gen- 
eral theoretical scenario of the vulcanization transition, 
but also they confirm the detailed analytical results for 
universal quantities, including the localized fraction ex- 
ponent and the distribution of scaled localization lengths. 



The ultimate origin of universality is not hard to un- 
derstand, despite the apparent intricacy of the order 
parameter associated with the amorphous solidification 
transition. As we saw in Sees. || and [II, there are two 
small emergent physical quantities, the fraction of local- 
ized particles and the characteristic inverse square local- 
ization length of localized particles. The smallness of the 
localized fraction validates the truncation of the expan- 
sion of the free energy in powers of the order parame- 
ter. The smallness of the characteristic inverse square 
localization length leads to a very simple dependence, 
X]q=o I'*"!^' n + 1 independent wave vectors of 

the replica theory, well beyond the permutation invari- 
ance demanded by symmetry considerations alone. As a 
result, near the transition, the amorphous solid state is 
characterizable in terms of a single universal function of 
a single variable, along with the localized fraction. 

Although throughout this Paper we have borne in mind 
the example of randomly crosslinked macromolecular sys- 
tems, the circle of ideas is by no means restricted to 
such systems. To encompass other systems possessing 
externally-induced quenched random constraints, such 
as networks formed by the permanent random covalent 
bonding of atoms or small molecules (e.g., silica gels), re- 
quires essentially no further conceptual ingredients (and 
only modest further technical ones) [p^ . 

One may also envisage applications to the glass transi- 
tion. Although it is generally presumed that externally- 
induced quenched random variables are not relevant for 
describing the glass transition, it is tempting to view 
the freezing-out of some degrees of freedom as the cru- 
cial consequence of the temperature-quench, with a form 
of quenched disorder thereby being developed sponta- 
neously. The approach presented in the present Paper 
becomes of relevance to the glass transition if one accepts 
this view of the temperature-quench, and thus models the 
nonequilibrium state of the quenched liquid by the equi- 
librium state of a system in which some fraction of cova- 
lent bonds has been rendered permanent (the deeper the 
quench the larger the fraction) ||l^. This strategy, viz., 
the approximating of pure systems by ones with "self- 
induced" quenched disorder, has also been invoked in 
very interesting work on the Bernasconi model for binary 
sequences of low autocorrelation |jl^. Interesting con- 
nections are also apparent with recent effective-potential 
approaches to glassy magnetic systems, in which one re- 
tains in the partition function only those configurations 
that lie close to the equilibrium state reached at the glass 
transition temperature pO| . 
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